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A quantum clock must satisfy two bàsic constraints. The first is a bound on the time resolution 
of the clock given by the difference between its maximum and minimum energy eigenvalues. The 
second follows from Holevo's bound on how much classical information can be encoded in a quantum 
system. We show that asymptotically, as the dimension of the Hilbert space of the clock tends to 
infinity, both constraints can be satisfied simultaneously. The experimental realization of such an 
optimal quantum clock using trapped ions is discussed. 
PACS numbers: 03.65.Bz, 42.50.Ar 



Recent technical advances in the làser cooling and 
trapping of ions suggest that coherent manipulat ions of 
trapped ions will be performed in the not too far fu- 
ture [Q . Apart from various important applications such 
as quantum information processing or improving high- 
precision spectroscopy these techniques also allow us to 
test fundamental concepts of quantum theory. In partic- 
ular, much deeper insight into the problem of quantum 
measurement can be obtained. 

In this Letter we study the problem of building an op- 
timal quantum clock from an ensemble of N ions. To be 
specific let us assume an ion trap with N two-level ions 
all in the ground state |\&) = |0) ® . . . <g> |0). This state 
is an eigenstate of the free Hamiltonian and thus cannot 
record time. Therefore the first step in building a clock is 
to bring the system to an initial state Cl which is not an 
energy eigenstate. For instance one can apply a Ramsey 
pulse whose shape and duration is chosen such that it 
puts all the ions in the product state 

Cl prod = p®...®p ; p=~(|0) + |l»((0|+(l|). (1) 

We shall also consider more general states, but shall al- 
ways take them to belong to the symmetric subspace of 
the N ions. The basis vectors of this space will be de- 
noted \m), m = 0, 1 . . . N. They are the completely sym- 
metrized states of N two-level ions with m ions in the 
ground state and (N — m) ions in the excited state. The 
states |m) have energy E m = m (this defines our unit of 
energy, setting h = 1 then defines our unit of time) . 

The reason we can restrict ourselves to the symmet- 
ric subspace is that we can map any clock state onto 
the symmetric subspace without affecting its dynamics. 
Indeed consider an initial state Cl — \%jj)(yj\ that does 
not belong to the symmetric subspace of the àtoms. Wc 
can decompose \ip) — ^2 m ^2 a c ma \m,a) where \m,a), 
a = 1, (^) denote a basis of the states with energy 
to. Consider the unitary operator U that maps the state 
\ip) onto the symmetric subspace without changing its 
energy: U^2 a c ma \m,a) = c m \m) where |to) is as beforc 
the symmetric state with energy to. Since U commutes 



with the Hamiltonian the performance of the clock based 
on \ip) is identical to the clock whose initial state is the 
symmetric state \tp S ym) = U\ip). 

After the preparation stage, the ions evolve in time ac- 
cording to the Hamiltonian evolution Cl(t) = U(t)ÚU' (t), 
U(t) = cxp{—ítH}. The task is to determine the elapsed 
time t by carrying out a measurement on the ions. Notc 
that because of the indeterminism of quantum mechanics 
it is impossible, given a single set of N two-level ions, to 
determine the elapsed time with certainty. The best we 
can do is to estimate the elapsed time based on the result 
of a measurement on the system . 

Making a good quantum clock requires a double opti- 
mization. First of all one can optimize the measurement. 
This aspect has been studied in detail in || where the 
best measuring strategy was derived. But one can also 
optimize the initial state Cl of the system. It is this sec- 
ond optimization that is studied in this Letter. 

Before turning to the problem of optimizing the ini- 
tial state Cl, it is instructive to review the fundamental 
limitations on the performance of quantum clocks. Let 
us first consider a simple classical clock that can then 
be generalized to the quantum case. Our classical clock 
consists of a set of n registers. Each register is either in 
the or the 1 state. Thus the classical clock consists of n 
bits, and can be in 2™ different states. The dynamics of 
the clock is as follows: the first register flips from to 1 
or from 1 to every 27r2~", the second register flips every 
27r2~™ +1 , etc... The last register flips every tt. This clock 
thus measures time modulo 2ir. Note that this clock has 
an inherent uncertainty since it cannot measure time with 
a precision better than 27r2~". Throughout this Letter 
the time uncertainty is defined as 



At 2 = [{t esUmate - í true )(mod 2ir)} 2 . (2) 

For the classical clock At c i ass = ir2~ n 

It is straightforward to replace the classical clock by 
a quantum version. The quantum clock consists of n 
two level systems (qubits). The first qubit has an en- 
ergy splitting between the levels of 2™ _1 so that it has 
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the same period as the corresponding classical register. 
The second qubit has an energy splitting of 2 n ~ 2 and so 
on up to the last qubit that has an energy splitting of 1. 
Considered together these n qubits constitute a quantum 
system with 2™ cqually spaccd energy levcls. 

The mapping between the Hilbert space of this ab- 
stract quantum clock and the symmetric subspace of N 
2-level ions is straightforward when N + 1 = 2 n . Indeed 
in this case the dimension and energy spectrum of both 
Hilbert spaces coincide. Note however that this is not 
a mapping between the qubits of the clock and the ions 
individually, but between energy eigenstates. This com- 
parison between classical and quantum clocks suggests 
that a quantum clock built out of N ions cannot be- 
have better then a classical clock built out of ln(iV + 1) 
registers. That this is indeed the case follows from two 
fundamental constraints: 

The first constraint is a bound on the time resolution 
of the clock that results from its energy spectrum. Indeed 
the time-energy uncertainty 



AtAE > 
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(3) 



where AE 2 = Tr(H 2 Ü) - [Tr(iíÒ)] 2 relates the uncer- 
tainty in the estimated time [defined by Eq. (||)] to the 
spread in energy of the clock. In the present case there 
is a state with a maximum energy uncertainty, namely 
\ip+) = + |0)) for which AE = N/2. Inserting in 

Eq (||) shows that for any clock built out of N ions 



(4) 



Note that one cannot attain equality in Eq. (||). Indeed 
the state evolves with a period 2tt/N, hence it nec- 
essarily has a large time uncertainty. Thus Eq. (||) is only 
a lower bound and maximising the energy uncertainty as 
in H is not necessarily an optimal procedure. 

The second fundamental constraint the clock must 
obey is a bound on the total information it can carry. 
Indeed Holevo || has shown that one cannot encode and 
subsequently retrieve reliably more than n bits of classi- 
cal information into n qubits. Letting a clock evolve for 
a given time interval t can be viewed as trying to encode 
information about the classical variable t into the state 
of the clock. Hence a measurement on our model clock 
(consisting of n — ln(7V + 1) qubits) cannot retrieve more 
than ln(N + 1) bits of information about t. 

Together these two bounds imply that the quantum 
clock cannot perform better than the corresponding clas- 
sical clock: it cannot carry more information and it can- 
not have better resolution. But is it possible, by making 
an optimal measurement and choosing in an astute man- 
ner the initial state of the clock, to make a quantum clock 
with similar performances to the classical one? Our main 
result is to show that this is indeed the case for clocks 
built out of symmetric states of N ions and to provide 
an algorithm for constructing such an optimal clock. 



The problem of constructing quantum clocks has been 
considered previously in fftj^. However the best clocks 
considered in these papers are based on the phase state 
l^o) described below. As we shall see for these clocks 
the time uncertainty is very large At ~ and is very 
far from reaching equality in Eq. (Q). Recently Vaid- 
man and Belkind ]9| considered the problem of clock for 
which equality holds in Eq. (^) . They showed that in the 
limit of large N the product states satisfy this condition. 
However for the product state the energy uncertainty is 
very large AE — yN /2 hence they also do not saturate 
Eq. (Q). Furthermore clocks based on product state also 
do not attain Holevo's bound. 

A similar approach to the one used here, namely opti- 
mizing both the initial state and the measurement on a 
system of N ions was considered in |ÏÏ| with the aim of 
using the ions as an improved frequeney Standard. This 
problem can be rephrased in the following way: one dis- 
poses of a classical but noisy clock which provides some 
a priori knowledge about the time t and one wants to 
improve the knowledge of t by using the N ions. On the 
other hand in the present Letter we suppose that there 
is no prior knowledge about t. The other difference with 
the present work is that our aim is to study the funda- 
mental structure of quantum mechanics. We therefore 
neglect the effect of noise during the preparation and 
measurement stages and decoherence during the evolu- 
tion. On the other hand taking these effects into account 
was central to |ï(| . 

In order to find how to build an optimal clock we must 
delve in detail into the functioning of this device. We 
first recali Holevo's results concerning the optimal mea- 
surement strategy ||. The measurement is described 
by a positive operator measurement (POVM), that is a 
set {O r }^ =1 of positive Hermitian operators such that 
J2 r O r = 1. To each outeome r of the measurement we 
associate an estimate t r of the time elapsed. The differ- 
ence between the estimated time t r and the true time t 
is measured by a cost function f(t r — t). Here we note 
that because of the periodicity of the clock, / has to be 
periòdic. We also take f(t) to be an even function. The 
task is to minimizc the mean value of the cost function 



2tt 



Tr[Ó r Ü(t)]f(t r -t)^. 

l'K 



(•5) 



To proceed we expand the cost function in Fourier series: 



f(t) = wq — Wk cos kt. 
fc=i 



(6) 



The essential hypothesis made by Holevo is positivity 
of the Fourier coefficients: Wk > 0, (A; = 1,2,...). He 
also supposes that the initial state Vl = \u>) = 

J2 m a m\m) is a pure state (and we make a phase conven- 
tion such that a m is real and positive). He then shows 
that 
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(7) 



k=l 



and equality is attaincd only if the measurement is of the 
form 



O r =p r \y r )(y r \ ; Pr >0 ; |* r ) -e lt ^|*o>, 

N 

(8) 



with the completeness relation O r = 1. 

Several remarks about this result are called for: 

1. Holevo supposed that the initial state is a pure 
state. If the initial state is mixed, Cl = 
^2iPi\tpi}{tpi\, then one finds that the correspond- 
ing cost is bounded by the average of the bounds 
Eq. (0). This shows that in building a good quan- 
tum clock one should always take the initial state 
to be pure. 

2. Holevo considcred covariant measurements in 
which the times t r takes the continuum of vàlues 
between and 2ir. But as shown in [jïïj the com- 
pleteness relation can also be satisfied by taking a 



discreet set of times í, 



2iL 
N+l ' ■> 



0,...,N. These 



"phase" states [jl2j form an orthonormal ba- 
sis of the Hilbert space, and this measurement is 
therefore a von Neumann measurement. It means 
that it is not necessary to use an ancilla to make 
the optimal measurement in this case. 

3. In Eq. (^) only the vàlues k = 0,...,N intervene 
because of the condition |m — m'\ = k. That is, 
only the first N + 1 Fourier coefficients of the cost 
function are meaningful. 

4. Because of the condition of positivity of Wk not all 
cost functions are covered by this result, but sev- 
eral important examples are: 4sin 2 | = 2(1 — cosí), 
|í(mod 2tt)\, |sin||, -<5[í(mod 2ir)}. The most no- 
table absent from this list is the quadratic devia- 
tion í 2 [as defined in Eq. (Q)] but it can be well 
approximated by the the first cost function since 
4sin 2 |~í 2 for |í| « n. 

We now turn back to the central problem of this Letter, 
namely how to optimize the initial state of the system so 
that the time estimate is as good as possible. This corre- 
sponds to minimizing the right hand side of Eq. (^) with 
respect to a rn . To this end let us define the matrix 



w 6 v 



1 N 

-Y 

2 ^ 

k=\ 



)• (9) 



We must then minimize / = a T Fa under the condition 
a T a = 1. Using a Lagrange multiplier we find the eigen- 
value equation (F — Xl)a = 0, and the task is therefore 
to find the smallest eigenvalue and eigenvector of F. 



Let us first consider the cost function / = 4sin 2 í/2. 
The advantage of this cost function is that the matrix F 
is particularly simple in this case. Furthermore for errors 
much smaller then 2tt, / and Aí 2 as defined in Eq. (||) 
coincide. Hence the first constraint on clock resolution 
Eq. (||) can be approximately replaced by / > 1/N 2 . 

If the initial state is the product state Eq. (Q) then 
the mean cost is given by the expression / = 2[1 — 

2~ Ar y^n 1 \ l( i v )( ÍLi)l which for large N decreases as 



J2i=0 Y ( i )( í+l! 

f ~ If the initial state is the phase state l 'í'o) j 

then direct calculation shows that / = j^p[- Thus for 

both states Aí ~ N~ i and one is very far from attaining 
equality in Eq. (||) . 

However neither of these two states is optimal. To 
find the optimum we note that the matrix F = 2S mm > — 
Smm'+i ~ ím+im' can be viewed as the discretized sec- 
ond derivative operator F ~ —d 2 /dx 2 with von Neu- 
mann boundary conditions. The lowest eigenvalue of F 
is therefore approximately X m in — 7r 2 /(iV + l) 2 and the 
corresponding eigenvector is 
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\- . 7^ + 1/2) , 

= > sin m) 

1^ N + l 1 1 

m 



(10) 



Thus in this case the cost decreases for large N as 



Opt — (AT + 1 )2 



corresponding to Aí 



opi 



{N+l)- 



There- 



fore, up to a factor of ir the optimal clock attains the 
bound Eq. (@). 

It is also interesting to consider the situation where the 
cost function is the delta function / = — <5[t(mod 2ir)]. In 



this case F„ 



2tt 



for all 



One checks that 



the phase state |\I/o) is the eigenvector of F with mini- 
mal eigenvalue A = — Note that one could also have 
taken a smeared delta function since only the first N + l 
terms intervene in Eq. (0). The smeared delta function 
is approximately zero everywhere except in an interval of 
about 1/N around zero where it is equal to N. Thus for 
this cost function one wants to maximize the frequcncy 
with which the estimated value of í is within about jj 
of the true value. But there is no extra cost if the esti- 
mated value is very far from the true one. It is for this 
reason that taking | ""í'o} as the initial state when the cost 
function is 4sin 2 í/2 is bad since making estimates that 
are wildly off is strongly penalized in that case. 

The mean value of a cost function gives only very par- 
tial information about the sensitivity of a clock. The 
full information is encoded in the probability P(í|í r ) = 
P(Ò r \t)P{t)/P{Ò r ) = P(Ó r |í)^±I that the truc time is 
í given that the read out of the measurement is t r . In 
the case of the optimal state \^ pt), one finds that 



P op t(t\t r ) — N- 



[1 + cos(iV+ Ï)T}(1 + cosT) 



1 - cos(T - 



N+l' 



1 - cos(T 



N+l' 



(11) 



3 



where M _ ^ 
tion (see Fig. |ÏJ) 



and T 



This distribu- 



has a central peak of width 3tt/(N + 1) 
and tails which decrease for 2ir >> \t — t r \ >> N^ 1 as 
P opt (t\t r )~N~ 3 \t-t r \~ 4 . 



For the phase state l^o) one finds 



P|* o) (í|í r ) = 



1 [1 -cos(AT+l)T] 
2tt(ÍV + 1) (1 -cosT) 



(12) 



This distribution has a slightly tighter central peak of 
width 2tï/(N + 1) but the tails of the distribution de- 
crease as iV -1 |í — t r \~ 2 . It is these slowly decreasing 
tails that give the main contribution to Aí ~ N~ 3 . 

For the product state Eq. (Q) the distribution 
Pprod(t\t r ) has a very wide central peak of width ~ -^=. 
In this case it is the wide central peak that gives rise to 
the large time uncertainty. 

Using these distributions it is possible to calculate the 
number of bits of information about time that is encoded 
in the outcomes of the measurement. This is given by the 
mutual information 

I=- J dtP(t)lnP(t) + ^P(i r ) J dtP{t\t r )\nP{t\t r ). 

(13) 

In all cases the integral in the second term is dominated 
by the central peak. Thus one finds that for the product 
state only i ln N bits of information about time are ob- 
tained, whereas for both the states |\&o) an( i V^opt) one 
obtains ln N bits, thereby saturating Holevo's bound. 

In summary we have seen that for different cost func- 
tions there are different optimal clocks. There is of course 
no cost function that is in an absolute sense better then 
anothcr, and the choice of a particular cost function de- 
pends on the physical context. Nevertheless the experi- 
mental realization of a quantum clock based on the state 
l'I'opt) seems particular ly desirable because it combines 
the attractive features that the a posteriori probability 
P(t\t r ) has a tight central peak and rapidly decreasing 
tails. 

Carrying out such an experiment with trapped ions 
presents two main difhculties. The first is the prepara- 
tion of the initial state Ú. Such coherent manipulation 
of trapped ions is one of the current experimental chal- 
lenges. A possibly important simplifying feature of the 
quantum clock is that since it is symmetric in the N 
ions one does not need to address each ion individually. 
The second problem is to measure the phase states 
Note that to perform such a measurement it is sufficicnt 
to be able to realize coherent manipulations of the ions 
and to measure the projectors \m){m\, i.e. to measure 
the number of excited ions. Precisely this type of the 
measurement has recently been realized ft3f . One may 
therefore hope that for systems with moderate number of 
trapped ions it will be feasible to make optimal quantum 
clocks in the not too far future. 
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FIG. 1. A plot of the a posteriori probability P(t\t r ) that 
the time was í given that the measurement yielded outeome 
t r = 7r for different initial states (N = 20). The dotted line 
corresponds to the product state Eq. (^), the sòlid line to the 
phase state \^o) and the dashed line to the optimal state 
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